SOS model partition function and the elliptic 
weight functions 

bo 

O ! S PakuliaktS, V Rubtsov«, A Silantyev^ 

t Laboratory of Theoretical Physics, JINR, 141980 Dubna, Moscow reg., Russia 
' t Department de Mathematiques, Universite d'Angers, 2 Bd. Lavoisier, 49045 

Angers, France 

" Institute of Theoretical and Experimental Physics, Moscow 117259, Russia 



< 

a: 

-I— > 



> 
in 

On 



O 

oo 
O 



X 



E-mail: pakuliak@theor.jinr.ru, Volodya.Roubtsov@univ-angers.fr, 
silant@tonton . univ-angers . f r 



Abstract. We generalize a recent observation [TJ that the partition function of the 6- 
vertex model with domain- wall boundary conditions can be obtained by computing the 
projections of the product of the total currents in the quantum affine algebra U q (sl2) 
in its current realization. A generalization is proved for the the elliptic current algebra 
[H[3]. The projections of the product of total currents are calculated explicitly and 
are represented as integral transforms of the product of the total currents. We prove 
that the kernel of this transform is proportional to the partition function of the SOS 
model with domain- wall boundary conditions. 
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1. Introduction 

The main aim of this paper is to apply the method of elliptic current projection to the 
computation of the universal elliptic weight functions. The projection of currents first 
appeared in the works of B. Enriquez and the second author [3], [H], as a method to 
construct a higher genus analog of the quantum groups in terms of Drinfeld currents [6]. 
The current (or "new" ) realization supplies a quantum affine algebra with a second co- 
product, the "Drinfeld co-product". The standard and Drinfeld co-products are related 
by a "twist" (see [4]). The quantum algebra is decomposed in two different ways a 
product of two Borel subalgebras. For each subalgebra, we can consider its intersection 
with these two Borel subalgebras and express it as their product. Thus we obtain for 
each subalgebra a pair of projection operators from it to each of these intersections. The 
above-mentioned twist is defined by a Hopf pairing of the subalgebras and the projection 
operators. See Section 4 where we recall an elliptic version of this construction. 

S. Khoroshkin and the first author have applied this method to a factorization 
of the universal i?-matrix [7] in quantum affine algebras, in order to obtain universal 
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weight functions [U [H] for arbitrary quantum affine algebras. The weight functions play 
a fundamental role in the theory of deformed Knizhnik-Zamolodchikov and Knizhnik- 
Zamolodchikov-Bernard equations. In particular, in the case of U q (gl n ), acting by the 
projection of Drinfeld currents onto the highest weight vectors of irreducible finite- 
dimensional representations, one obtains exactly the (trigonometric) weight functions 
or off-shell Bethe vectors. In the canonical nested Bethe Ansatz, these objects are 
defined implicitly by recursive relations. Calculations of the projections are an effective 
way to determine the hierarchical relations of the nested Bethe Ansatz. 

It was observed in [1] that the projections for the algebra U q (sl 2 ) can be represented 
as integral transforms and that the kernels of these transforms are proportional to the 
partition function of the finite 6-vertex model with domain-wall boundary conditions 
(DWBC) pp. We prove that the elliptic projections described in [2] make it possible 
to derive the partition function for elliptic models. We show that the calculation of 
the projections in the current elliptic algebra [21 [3] yields the partition function of the 
Solid-On-Solid (SOS) model with domain-wall boundary conditions. 

The partition function for the finite 6-vertex model with domain wall boundary 
conditions was obtained by Izergin [9], who derived recursion relations for the partition 
function and solved them in determinant form. The kernels of the projections satisfy 
the same recursion relations and provide another formula for the partition function. 

The problem of generalizing Izergin's determinant formula to the elliptic case has 
been extensively discussed in the last two decades. One can prove that the statistical sum 
of the SOS model with DWBC cannot be represented in the form of a single determinant. 
While this paper was in preparation, H. Rosengren [10] showed that this statistical sum 
for an n x n lattice can be written as a sum of 2 n determinants, thus generalizing 
Izergin's determinant formula. His approach relates to some dynamical generalization 
of the method of Alternating-Sign Matrices and follows the famous combinatorial proof 
of Kuperberg [16] . 

We expect that the projection method gives a universal form for the elliptic weight 
function [TJ] as it does for the quantum affine algebras [12]. When this universal weight 
function is represented as an integral transform of the product of the elliptic currents, 
we show that the kernel of this transform gives an expression of the partition function 
for the SOS model. On the one hand we generalize Izergin's recurrent relations and 
on the other hand we generalize to the elliptic case the method proposed in [1] for 
calculating the projections. We check that the kernel extracted from the universal 
weight function and multiplied by a certain factor satisfies the recursion relations that 
have been obtained, which uniquely define the partition function for the SOS model 
with DWBC. Our formula given by the projection method coincides with Rosengren's. 

An interesting open problem which deserves more a extensive study is the relation 
of the projection method with the elliptic Sklyanin-Odesskii-Feigin algebras. It was 
observed in the pioneering paper [3] that half of the elliptic current generators satisfy the 
commutation relations of the VT-elliptic algebras of Feigin. Another intriguing relation 
was observed in [17] : there existe a certain subalgebra in the "A-generalization" of the 
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Sklyanin algebra such that its generators obey the Felder's .R-matrix quadratic relations 
given in [IB] . The latter paper gives also a description of the elliptic Bethe eigenvectors 
(the elliptic weight functions). 

This is a strong indication that the projection method should be considered and 
interpreted in the framework of the (generalized) Sklyanin-Odesskii-Feigin algebras. We 
hope to discuss this problem elsewhere. 

The main results of this paper were reported at the 7-th International Workshop on 
"Supersymmetry and Quantum Symmetry" in JINR, Dubna (Russia), July 30 - August 
4, 2007. 

The paper is organized as follows. In section [2] we briefly review the finite 6- vertex 
model with DWBC, and we present the formulae for the partition function: Izergin's 
determinant formula and the formula obtained by the projection method. Section [3] is 
devoted to the SOS model with DWBC. We briefly introduce the model and pose the 
problem of how to calculate the partition function of this model. We derive analytical 
properties of the partition function that allow us to reconstruct the partition function 
exactly. In section H] we introduce the projections in terms of the currents for the 
elliptic algebra, following [2]. We generalize the method proposed in [1] to this case in 
order to obtain the integral representation of the projections of products of currents. 
Then, using a Hopf pairing, we extract the kernel and show that it satisfies all the 
necessary analytical properties of the partition function of the SOS model with DWBC. 
In Section [5j we investigate the trigonometric degeneration of the elliptic model and of 
the partition function with DWBC. We arrive at the 6-vertex model case in two steps. 
The model obtained after the first step is a trigonometric SOS model. Then we show 
that the degeneration of the expression derived in Section H] coincides with the known 
expression for the 6-vertex model partition function with DWBC. An appendix contains 
the necessary information on the properties of elliptic polynomials. 

2. Partition function of the finite 6-vertex model 

Let us consider a statistical system on a square n x n lattice, where the columns and 
rows are numbered from 1 to n from right to left and from bottom to top, respectively. 
This is a 6-vertex model where the vertices on the lattice are associated with Boltzmann 
weights which depend on the configuration of the arrows around a given vertex. The six 
possible configurations are shown in Fig. [TJ The weights are functions of two spectral 
parameters z, w and ananisotropy parameter q: 

a(z,w) = qz — q~ 1 w, b(z,w) = z — w, 

c(z, w) = (q — q~ l )z, c(z, w) = (q — q~ )w. 

Let us associate the sign '+' to the arrows directed upward and to the left, while 
the sign ' — 'is associated to to the arrows directed downward and to the right as shown 
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Figure 1. Graphical presentation of the Boltzmann weights. 



in Fig. [TJ The Boltzmann weights ([T 

( a(z,w) 

R(z, w) = 
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are gathered in the matrix 
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b(z, w) c(z,w 

c(z,w) b(z,w 




a(z,w) J 



(2) 



\«/3 



acting in the space C 2 ® C 2 with the basis e a ® e^, a,/3 = ±. The entry R(z,w)^ s , 
a, [3, 7, 5 — ± coincides with the Boltzmann weight corresponding to Fig. [2j 
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Figure 2. The Boltzmann weight R(z,w)"g. 

Different repartitions of the arrows on the edges form different configurations 
{C}. A Boltzmann weight of the lattice is the product of the Boltzmann weights 
at each vertex. We define the partition function of the model as the sum of the 
Boltzmann weights of the lattice over all possible configurations, subject to some 
boundary conditions: 



2(M,M) 



{C} i,j=l 



(3) 



Here a^, fy, jij, 5ij are the signs corresponding to the arrows around the (i, j)-th 
vertex. We consider an inhomogeneous model where the Boltzmann weights depend on 
the column by the variable Zi and on the row by the variable Wj (see Fig. [3]). 

We choose the so-called domain-wall boundary conditions (DWBC) that fix the 
boundary arrows (signs) as shown in Fig. [3j In other words, the arrows enter on the left 
and right boundaries and leave on the lower and upper ones. 
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Figure 3. Inhomogeneous lattice with domain wall boundary conditions. 



In [9], A. G. Izergin found a determinant representation for the partition function 
of the lattice with DWBC, 



z ({ z }A w }) = (?-?' 



-l\n 



n 

m=l 



W m X 



X- 



I J (Zi - Wj)(qzi - q 1 w j ) 
»J=1 



det 



1 



[Zi - Wj)(qZi - q hUj) 



(4) 



t,i=l,...,n 



n (* 

n>i>j>l 

Izergin's idea was to prove a symmetry of the polynomial ([3]), and then use it to find 
recursion relations for the quantity Z({z}, {w}) and to observe that these recursion 
relations allow the reconstruction of Z({z}, {w}) in a unique way and that the same 
recursion relations are valid for the determinant formula (j3J). 

On the other hand it was observed that the kernel of the projection of n currents is 
a polynomial of the same degree, and satisfies the same recursion relations pQ. It means 
that this kernel coincides with the partition function for the n x n lattice. Moreover, 
the theory of projections gives another expression for the partition function: 



n 



z {{ z )A w }) = (q-q 



-l\n 



W m } f 

m=l n>i>j>l 



n 



q Wi — qwj 



x 



Wi — w 



x 



e n 

crgS„ l<i<j<n 



-1 , 



qWg(i) ~ q ^Wgjj) 

q^w a(i) - qw u{j) n ; l j 



] [ (qzi-q 1 w u{k ) j [ (Zi-w a ^), 



(5) 



l<i<k<n 



q~ 1 ™<T(i)-qu>a(j) 

product if both conditions i < j and a(i) > a(j) are satisfied simultaneously. 



where S n is the group of permutations. Here the factor H "^ w H _ UJaU) appears in the 



3. Partition function for the SOS model 

3.1. Description of the SOS model 

The SOS model is a face model. We introduce it in terms of heights as usual, but then 
we represent it in the R- matrix formalism as in [13]. This language is more convenient to 



SOS model partition function and the elliptic weight functions 



6 



n 



3= 2 
1 


i= n ... 2 10 
Figure 4. The numbering of faces. 



generalize the results reviewed in Section [2] and to prove the symmetry of the partition 
function. 

Consider a square nxn lattice with vertices enumerated by the index i — 1, . . . , n as 
in the previous case. It has (n+1) x (n+1) faces enumerated by pairs i, j = 0, . . . , n 
(see Fig. H]). To each face we assign a complex number called its height in such a way 
that the differences of the heights corresponding to the neighboring faces are ±1. Let 
us denote by dij the height corresponding to the face placed to the upper left of 
the vertex (i,j). Then the last condition can be written in the form [dy — di-ijl = 1, 



0, . . . , n, and \dij 



d 



ij-l\ 



1 for i = 0, . . . , n, j 



1, . . . , n. Each 



for % = 1, . . . , n, j 

distribution of heights d^ — 0, . . . , n) subject to these conditions and to boundary 
conditions defines a configuration of the model. It means that the partition function of 
this model can be represented in the form 



II W ij(di,j-udi 

—l,j—l: di—l,jj dij): 

C i,j=l 



(6) 



where Wij(dij-i,di-i t j-.i,di-.i ) j,dij) is the Boltzmann weight of the (i,j)-th vertex 
depending on the configuration by means of the heights of the neighboring faces as 
follows \U 



Wij(d + l,d+ 2, d+ l,d) 
Wij(d-l,d-2,d-l,d) 



a[Uj 



b(ui 




M) 


b(ui 


-vf, 


hd) 


c(Ui 




hd) 


c(Ui 


-Vj; 


Hd) 



6(ui 


- vj + h), 


e(ui 


- Vj + h), 


0{ui 


- Vj)6{hd + h) 




6{hd) 


B{ui 


- V j)6{hd - h) 




6{hd) 


0(ui 


- Vj + hd)6(h) 




6{hd) 




- Vj - hd)9(h) 


6{-hd) 



(7) 



Wij(d+ l,d,d + l,d) 

As in the 6- vertex case the variables Ui, Vj are attached to the i-th vertical and j-th 
horizontal lines respectively, h~ is a nonzero anisotropy parameter 0. The weights are 

| In the elliptic case, we use additive variables Ui, Vj and an additive anisotropy parameter h instead 



of the multiplicative variables Zi 



and the multiplicative parameter q = e" 
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Figure 5. The Boltzmann weights for the SOS model. 



expressed by means of the ordinary odd theta-function defined by the conditions 



6(u + 1) = -9(u), 6(u + t) 
Let us introduce the notations 



-2iriu— irir 



9(u), 0'(O) = I. 



(8) 



01 ij 



di 



di 



i—±j ^iji Ptj di— ij— i di_ij, di_\j_\ d{j_i, 5ij dij—% d^j. (9) 

The differences ()9]) take the values ±1 and we attach them to the corresponding edges 
as in Fig. [2j jij+x = a.%j is the sign attached to the vertical edge connecting the 
th vertex to the + l)-st one, Aj+i = $ij is f ne s ig n attached to the horizontal 
edge connecting the (i,j)-th vertex to the (z + 1, j)-th one. The configuration can be 
considered as a distribution of these signs on the internal edges subject to the conditions 
ocij + Pij = jij + Sij, i, j = 1, . . . , n. In terms of signs on the external edges the DWBC 
are the same as shown in Fig. [3j Additionally, we have to fix one of the boundary 
heights, for example, d nn . 

The Boltzmann weights (j7|) can be represented as the entries of a dynamical elliptic 
i?-matrix [131: 



Wij(di 7 j-x, dj-ij-i, d- 
R(u; A) 

V 



(10) 



R{ui 

( a(u) \ 

b(u;X) c{u;X) 
c(u;X) b{u;X) 
a(u) ) 

Let T"™ (ui, {v}, Xi) be the column transfer matris. It is a matrix-valued function of 
Ui, all spectral parameters Vj, j = 1, . . . , n and the parameters Xi related to the heights: 

X^(«i,W,A i )fc^= (11) 



R^' n \ui - V n , A m )^ (n+1 ' n " 1} (w t " «n-U \n-l) ■ ■ ■ R^Hui - Vl ; Ail)) 



!Pin ■ ••Pi'. 
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R(n+^)( Ui _ Vn] K in )R^+^){ Ui - « n _ i; A i>n _i) • • • R^\ Ul - v r , Aa) 



1 



where Ay = My = A* + K S ih \ = hd in = A + h J2 a in, Kj = \ + h H^. 

l=j + l l=i+l 1=3+1 

The matrix acts in the Z-th two-dimensional space V\ = C 2 as diag(l, —1) and the 
i?-matrix R( a ' b ^ acts nontrivially in the tensor product V a (g> H- The superscript n + 1 in 
the /^-matrices is regarded as belonging to an auxiliary space V n+ i = C 2 . The partition 
function (jBJ) corresponding to DWBC (a>i n = +1, (3u = — 1, 70 = —1, 5 n i = +1, 
i — 1, . . . , n) can be represented by means of the column transfer matrices: 

Zl+({u}, {v}; A) = (T±(u u {v}, AO • ■ ■ T±(u n , {v}, A*))"""", (12) 

where Aj = A + fr(ri — z) . Similarly one can define the row transfer matrix. 

3.2. Analytical properties of the partition function 

Here we describe the analytical properties of the SOS model partition function which 
are analogous to those used by A.G. Izergin in order to recover the partition function 
of the 6-vertex model. These properties uniquely define this partition function. 

Proposition 1 The partition function with DWBC Z_+({u}, {v}; A) is a symmetric 
function in both sets of variables Ui and Vj. 

The proof is based on the dynamical Yang-Baxter equation (DYBE) for the R- 
matrix [13] 

R^\h - t 2 ; A)i? (13) (ti - t 3 ; A + hH^)R^\t 2 - t 3 ; A) = 

= R^\t 2 -t 3 ;X + hH^R^ih - t 8 ; A)i? (12) (ti - t 2 ; A + KH&). 

In order to prove the symmetry of the partition function Zt+({u}, {v}; A) under the 
permutation Vj Vj-i, we rewrite the DYBE in the form 

R( n+1 >i\ Ui - v 3 - A^R^-^im - A,, + hH&) 

x RM-Qfa - v^i, Aij) = /,"••' l: (r.; - Vj -n Ay + hH^) (13) 

x J&+ x J- x \ui - Ay)# (n+1,j) (ui - uj; Aij + HH^-V). 

Multiplying the z-th column matrix f fTTT) by R^^~ l \vj — Uj-ijAy) to the right and 
moving it to the left using (TT31 . the relation [f/i + H 2 ,R(u, A)] = and the equality 
Ay + ^«j n = Aj_ij, we obtain 

T^( Wi , {v}, \;)ir ! - ' iV ; - Ay) = ' : ir ; - Vj.!] Aj_ij) 

xp^- 1 )^™^, <- A,) ■ • ■ t;7( M(1) <- ^J, K)V^- X \ (14) 

where "P G End(C 2 ®C 2 ) is a permutation matrix: P(ei<8>e2) = e2<S>ei for all ex, e 2 G C 2 
and the notation {u,- Uj-i} means that the set of parameters {v} with Vj^i and i>j are 
interchanged. Multiplying the product of the column matrix by R^^~ l \vj — fj-i; A n j) 
to the right and moving it to the left using (Tbil) one yields 

T«-K {„}, AO • • • K, {«}, \„)l! ri l \ Vj - A ni ) = - A 0J .) 

xV^T^im, { Vj <- Ai) • • • t^;;K, - ^J, K)v^~ x \ (15) 
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n n 

where A j = A + cti n + h~ Yl H^ l \ A n j = X n = X. Finally, comparing the matrix 
i=l i=j+i 

element (■)+'" + of both sides of ffT5l) . taking into account Formula f|T2|) and the identities 

R(u,X)~ s ~ = a(u)5-5s, R(u, A)f + = a{u)^, Pj- = 6;5j, V ++ = 

(where 8® is the Kronecker symbol) and substituting a in = +1, jn = —1, one derives 

Z±+{{u}, {«}; X) = Zt-({u}, { Vj «-> v^}; A). (16) 

Similarly, using the row transfer matrix one can obtain the following equality from the 
DYBE: 

Z±+({u}, {v}; A) = Z±~({ Uj <- {«}; A). (17) 

The partition function with DWBC satisfies relations ( fTB"l) . (1171) for each j = 1, . . . ,n, 
which is sufficient to establish the symmetry under an arbitrary permutation. □ 

Proposition 2 The partition function with DWBC is an elliptic polynomial [jj of 
degree n with character x in each variable Ui, where 

n 

X(l) = (-l) n , X(r) = (-1)" exp (2tu(A + £ Vj j) . (18) 

i=i 

Due to the symmetry with respect to the variables {u} it is sufficient to prove 
the proposition for the variable u n . To represent explicitly the dependence of 
{{u}, {v }; A) on u n , we consider all the possibilities for the states of the edges 
attached to the vertices located in the n-th column. First, consider the (n, n)-th 
vertex. Due to the boundary conditions a nn = 5 nn = +1 and to the condition 
ctnn + Pnn = Inn + &nn we have two possibilities: either (3 nn = ^ nn = —1 or 
Pnn = Inn = +1- In the first case, one has a unique possibility for the rest of the 
n-th column: j n j = —1, (3 n j = +1, j = 1, ...,n — 1; in the second case, there 
are two possibilities for the (n,n — l)-st vertex: either (3 n>n ^\ = 7„ in _i = —1 or 
Pn,n-i = 7n,n-i = +1; etc. Finally the partition function is represented in the from 

n n 

z -+(i u }i M; A) = Yl a{u n - Vj) c(u n - v k ; X + (n - k)H) 
k=i j=k+i 

jfc-i 

x o(w n - vj; A + (n - j)h)g k {u n -i, . . . , «i, {f }; A), 
i=i 

where gk(u n -i, ■ ■ ■ > u i> {^j; A) are functions which do not depend on w n . Each term of 
this sum is an elliptic polynomial of degree n with the same character ( ITBl) in the variable 
u„. □ 

Remark 1 Similarly, one can prove that the function Z_+({u}, {v}; A) zs an elliptic 
polynomial of degree n with character x in each variable v i; where xiX) = ( — l) n ; 
x{t) = (-l) n e 27ri(_A+ ^=i" l) . 

§ The definition of elliptic polynomials and their properties are given in |Appendix A| 
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Proposition 3 The n-th partition function with DWBC ( f7l|j with the condition 
u n = v n — h can be expressed through the (n — l)-st partition function: 

Zt+{u n = v n - h,u n -i, . . . ,Mi;f n ,w n _i, ...,vx',X)= (19) 

6(\ + nh)6(h) ^fj ( <\r 7 +- ( u 

= q(\ + ( n - l)h) -1-1 V ' " ~ Vm ~ h > e ( Um ~ Vn > ) z -+( u n-i, ■ ■ • ,«i;u n -i, • • • A). 

Considering the n-th column and the n-th row and taking into account that 
a(u n — v n )\ Un=Vn _ h = a(—h) = we conclude that the unique possibility for a non- 
trivial contribution is: (3 nn = 7„ n = -1, 7„j = -1, (3 nj = +1, j = l,...,n — 1, 
(3in = —1, 7j n = +1, i = 1, ... ,n — 1. The last formulae impose the same DWBC for 
the (n — 1) X (n — 1) sublattice: 5 n -i,j = Ay = +1, J = 1,...,?^ — 1, a^n-i = 7m = +1, 
i = 1, ... ,n — 1, d n _i in _i = d nn . Thus the substitution u n — v n — H to the partition 
function for the whole lattice yiefds 

Z±+{u n = v n - H, M n _i, ...,Ui, v n , W n _i, ...,Vi,X) 
n—1 n~l 

= c(-H; A) 6(f n - Vj - h; A + (n - j)h) 6(u; - A + (n - 

xZ±+(u n -i, . . . ,mi;w„_i, . . . ,vt; A). (20) 

Using the explicit expressions ([7]) for the Boltzmann weights, one can rewrite the last 
formula in the form (fT9|) . □ 

Remark 2 From Formula UZDi) we see that the following transformation of the R-matrix 

b(u,v;\)^pb(u,v;\), b(u, v ; A) — ► p^ 1 b(u, v ; A) (21) 

does not change the recursion relation / TJPj) . where p is a non-zero constant which does 
not depend on u, v and A. 

Lemma 1 // the set of functions {Z^ l \u n , . . . , u\\ v n , . . . , V\\ A)} n >i satisfies the 
conditions of Propositions CI and the initial condition 

7(1), U -/ \ OjU! - V! - \)6(h) 

Z y >{u\\ v\\ A) = c[u\ - vx) = — — (22) 

then 

Z±~(u n , ...,ux;v n ,.. .,vx; A) = Z {n \u n , . . . , u X ] v n , . . .,v x \ A). (23) 

Due to (1221) . this lemma can be proved by induction on n. Let the equality (1231) 
be valid for n — 1. Consider the functions Z_+(u n , . . . , u X ', v n , . . . , Vx', A) and 

Z (n) 

(u n , . . . , ux, v n , . . . , Vx; A) as functions of u n . Both are elliptic polynomials of degree 
n with character ffTSl) . They have the same value at the point u n = v n — h~, and due to 
the symmetry of these functions with respect to the parameters {u,}™ =1 they coincide 
at all points u n = Vj — h, j = 1, . . . , n. It follows from Lemma [2] (see Appendix A ) that 
these functions are identical. □ 

Remark 3 As we can see from the proof of Lemma [21 it is sufficient to establish the 
symmetry only with respect to the variables Vj. 
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Remark 4 The transformation < f^7]j of the R-matrix does not change the partition 
function with DWBC. 

4. Elliptic projections of currents 

Let /Co = [[«]] be the completed space of complex-valued meromorphic functions 

defined in the neighborhood of the origin which have only simple poles at this point. 
Let {e*} and {cj} be dual bases in /C such that § ^-^(u) €j(u) = <5!-. 

4-1. Current description of the elliptic algebra 

Let A be a Hopf algebra generated by elements h[s], e[s], f[s], s G /C , subject to the 
linear relations 

x[a\Si + a 2 S2] = aix[si] + a 2 £[s 2 ], ai,a 2 eC, s±, s 2 6 /Co, 
where x G {h, e, /}. The commutation relations will be written in terms of the currents, 



The currents e(u) and f(u) are called the total currents. They are defined in terms of 
dual bases of )C and their definition does not depend on the choice of these dual bases 
(see [HI [15] ) . The currents h + {u) and h~(u) are called the Cartan currents and they are 
defined in terms of the special basis 



h + (u) = J2h[^°hAn), h-(u) = -5>fe;o]e ii0 («) 






(24) 






[K + (u),K-(v)] = 0, 



-l 




6{u -v + h) 



6{u — v — h)e{u)e{v) 



-i 



6(u-v- h) 

e(u-v + h) J[ * vh 

6{u — v + h)e(v)e(u), 



(25) 



6(u - v + h)f(u)f(v) = 9(u-v- h)f(v)f(u), 



(26) 
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where K + {u) = exp ( ^^23 ~~ h + (u) ) , K~{u) = exp(hh~(u)) and 5(it, ?;) = is 

a delta-function [jjj for /Co- The algebra A is a non-central version of the algebra A{t) 
introduced in [3]. This algebra is equipped with the co-product and co-unit: 

AK ± (u) = ® K^u), 

Ae(«) = e(«) (g) 1 + #~(tt) ® e(w), 

Af(u) = f(u)®K + (u) + l®f(u), 

e(K ± (u)) = l, e(e(u)) = 0, e(/(u)) = 0. 

Let Af and Ae be the subalgebras of A generated by the generators /i[ej ; o], f[s], 
and /i[e ?;0 ], e[s], respectively, s G /Co- The subalgebra Af is described by the currents 
K + (u), f(u), and the subalgebra Ae by K~{u), e(u). We introduce the notation H + 
for the subalgebra of A generated by ^[e»;o] > As stated in [2], the bialgebras (Af, A op ) 
and (.Ae, A) are dual with respect to the Hopf pairing (-,•}: A f x Ae — > C defined in 
terms of currents as follows: 

0(u — v — K) 
0{u-v + h) 

These formulae uniquely define a Hopf pairing on Af x Ae- In particular, one can 
derive the following formula 



(/(u),e(«)) = h-Hfav), (K+(u),K-(v)) = lf u _ V K l (27) 



(f(t n )---f(t 1 ),e(v n )---e(v l ) 



<t(;)><t(;') 

^.S. Projections of currents 

We define the projections as linear maps acting on the subalgebra Af- Dual projections, 
which we do not consider here, act in the subalgebra Ae- We define the projections in 
terms of the half-currents fx( u ) an d f\(u), defined below. These are usually defined 
as parts of the sum (12^1) (with the corresponding sign) such that f(u) = f\(u) — f\(u)- 
Here A is the parameter for the decomposition of the total current into the difference 
of half- currents. Elliptic half-currents are investigated in details on the classical level 
in [15]. We will introduce the half-currents by their representations by means of integral 
transforms of the total current f(u): 

, N f dv d(u-v - X) I dv 6{u-v - X) 

K {U)= f 2tH 9(u - v)9(-X) fAu)= f 2^ 6(u-v)6(-X) m > (29) 

|u|<|it| |d|>|« 

where A ^ T = Z + Zr. The half-current f\(u) is called positive and f\(u) is called 
negative. 

|| One can find more details about distributions acting on ICq and their significance in the theory of 
current algebras in our previous paper [15] . 
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The corresponding positive and negative projections are also parameterized by A 
and they are defined on the half-currents as follows: 

Wa + W) = Wa») = 0, (30) 

(«)) = 0. Wa») = /*(«)■ (31) 
Let us first define the projections in the subalgebra Af generated by the 
currents f(u). As a linear space this subalgebra is spanned by the products 
/( w n)/(Mn-i)---/(«i), n = 0,1,2,.... It means that any element of Af can be 
represented as a sum (maybe infinite) of integrals 

~7k W~~f{ u n) ■ ■ ■ f(ui) S n {u n ) ■ ■ • Si(ui), S n ,...,S 1 E Kq. 

(2m) n 

It follows from the PBW theorem proved in [2] that any element of Af can also be 
represented as the sum of the integrals 

• • • fx + 2mH( U m+l)ft+2(m~l)h( U m) ' ' ' /a + («i) S »W ' ' " 

s n ■ ■ - Si G /Co, < m < n. Therefore, it is sufficient to define the projections on these 
products of half- currents: 

P+(x~x + ) = e(x-)x + , P^(y-y + ) = y-e(y + ), (32) 

where 

X = f\+2(n-l)h( U n) ' ' ' f\+2mh( U m+l), V = f\ ( u n) ' ' ' /x_ 2 ( n _ m _l)ft( M m+l) ' 

X+ = f\+2(m-l)fr[ U ™) ' ' ' f\( U ^)i V + = f\-2(n~m)h( U m) ' ' ' f X-2(n-l)h( Ul ) ' 

The product of zero number of currents is identified with 1 and in this case: e(l) = 1. 
The counit e of a nonzero number of half-currents is always zero. So, this definition 
generalizes Formulae (1301) and ( 13TT) . We complete the definition of the projections on 
the subalgebra Af — Af ■ H + by the formulae 

P+(at + ) = P A + (a)t + , P A "(at+) = P^a)e(t+), 

where a G .4./, t + G H + . 

4-3. The projections and the universal elliptic weight function 
Consider the expressions of the form 

K-(n-l)n(f( U n)f(u n -l) ' ' ' fMfM), (33) 

where the parameter A — (n — is chosen for symmetry reasons. Let us begin with the 
case n — 1. Formula ( 1301) implies that in this case the projection is equal to the positive 
half-current, which can be represented as an integral transform of the total current: 



27T10(« 1 - vi)0(-A)' 
l«ll>[«ll 

The integral ^ without limits means a formal integral - a continuous extension of the integral over 
the unit circle. 



SOS model partition function and the elliptic weight functions 14 

The kernel of this transform gives the initial condition for the partition function with a 
factor: 

ZU( Ul ; Vl ; A) = e{h)e{ Ul - vi) ft d ( {Ul ~ V !7 ( X \ y (34) 

The projections (1331) can be calculated by generalizing the method proposed in [1] 
for the algebra U q (si2). The method uses a recursion over n. Let us first present the 
last total current in fl33l) as the difference of half- currents: 

Pt { n-i ) nU{u n )---f{u 2 )f{u 1 ))= (35) 

= P \-(n-3)h(f( U n) ■ ■ ■ f(u 2 ))ft(n-l)n( U l) ~ P X-(n-l)h(f 'K) " " " f (u 2 ) fx- [n -l)h( U l)) ■ 

In the first term we move out the positive half-current from the projection and, therefore, 
calculation of this term reduces to the computation of the (n — l)-st projection. In the 
second term in (13"oT) we move the negative half-current to the left step by step using the 
following commutation relation [2] 

ft \*-t \ 9{p-ui-K) 0(v - ui + A + H) 



where 



FM = -Q^^UU^ftiv) - fx+2h( v )fx( v ))- 



At each step we obtain an additional term containing F\(u) and at the last step 
the negative half-current is annihilated by the projection: 

n 

Px-(n- m (f(Un) ■ ■ ■ f(u 2 )fx_ in _ m (Ul)) ^Q^Xj, (36) 

3=2 



where 



n , x 9(uj -u + \-(n-2j + 2)K) i4 6{u k -u-H) 



9(uj — u + H) 0{u k — u + ft) ' 

X 3 = P X-(n-l)h(f( U n) ■ ■ ■ /(u i+ i)F A _ (n _j y+3 ) R (u i )/(u i _ 1 ) • • • f(u 2 )). 

Setting u\ = U{ in (1361) . we can substitute the negative half-current for the positive 
one using the commutation relation for the total currents f(u) and the equality 
f(u)f(u) = 0. Moving out the positive half-current to the left one obtains a linear 
system of equations for Xi, i = 2, . . . , n: 

n 

Px- { n-z)nU{u n ) ■ ■ ■ f(u 2 ))ft {n _ 1)h (u t ) = QjMXj. (37) 
Multiplying each equation fl37j) by 



6(Uj -11 + X) -p-r 8(ll k -Uj + tl) -p-r 8(u k - U) 

0(A) 11 9{u k -u + h) 11 9{u k - Ui ) ' 
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summing over % — 2, . . . , n and using the interpolation formula (see Appendix A) 

qau) = g Qi("i) m n n ^5 (^) 

yields 

P A + -(n-3),(/K) • ■ " /(«,)) g — ^ A fl( Mfc - « + fi) S ^3^) A-(n-l) ft (^) 

71 

= ^g i («)X J , (39) 
Comparing ( |39l) with ( |36i) . we conclude that 

^A + -(„-l )?i (/K) ■ ■ • /M/^^l)) = P Un-^U{Un) ■ ■ ■ f(u 2 )) 
(Ui - U\ + A) JT 0(li A - Wj + ^) TT 0(«fc - «l) 



X 



E PV"z - "1 -t- A) TT v\,Uk - <H -T ii) y-r oyu k - up + 
9{\) 11 - «! + n) 11 0(«* - ui) ^-(-d^- ^ 4U ^ 



i=2 v 7 fc=2 v x ' k=2 



Finally, returning to Formula (1351) we derive the following expression for the projection 

P X-(n-l)h(f( U n) ■ ■ ■ /(«2)/(«l)) = ^A + -(n-3)ft(/K) ' ' ' f i^)) ft(n-l)h( U ^ U ^ ■ ■ ■ > « 2 ), (41) 

where we introduce the linear combination of the currents: 

,+ / x _ , + ( \ _ ®( Ui - w. m + A + (m - l)fi) 

^A-(n-2m+l)av Um ' " " " ' U m+l) ~ J X-(n-2m+l)H\ Um ) 2-U 6( X + (m — 1)J%) 

i=m+l ^ ^ ' ' 

A flfofc -Uj + h) A g(Ufc - Urn) + . x f x 

ill - M ™ + Jl - «*) h ^- 2m+ ^ [Ul) - {42) 

Continuing this computation by induction we obtain an expression for the projections 
in terms of the half-currents (T42l): 



P X-(n-l)h(f( U n) ■ ■ ■ fMfM) = Y[ fx-(n-2m+l)h( U m; ««,•••, Vl)' ( 43 ) 

n>m>l 

Using the addition formula 

n n n 

n g Xi ( Ui - v) = ^2 n °^ («* - ^)^o («< - «), (44) 
i=i i=i j=i 

n 

where G\(u — v ) = jS^, Ao = X] one can represent the half-currents (1421) as 

j=i 

integral transforms of the total current: 

f\-{n-2m+l)hi Um 'i Un i • • • ' U m+l) = 

n O(u k - u m ) J" dv 6(u rn — v — \ — {m — l)h) A 9{u k - v + h) 
. M e( Uk -u m + h) f 27rie(u m -v)e(-x-(m-i)h) 11 0(« fc -«) ;w " 

fe=m+l i i. i i fc=m+l 
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Replacing each combination of the half-currents (142 p in (T4"5|) by their integral form we 
obtain 

D + (t( \ ft \tr w TT 9(u k -u m ) I dv n ---dvi 
P A + W/W-/W^))= II 9(u k -u m + ft) f ~12^- 

n>k>m>l \ui\>\vj\ 

n d{u k -v m + ft) 0{u m — f m — A — (m — l)ft) 

n>k>m>l m=l x ' v 

Formulae (|43p and (|45p yield expressions for the universal elliptic weight functions 
in terms of the current generators of the algebra A. 

4-4- Universal weight function and SOS model partition function 

To extract the kernel from the expression (T4"5l) and derive a formula for the partition 
function we use the Hopf pairing (1271) . Let us calculate the following expression 
generalizing (JMj) : 

ryfn), w TT at \ TT 0(u k - u m + ft)6(v k - v m - ft) 

Z< >(„„, Ml ; , i; A) = n »(«. " «,) I! e { u k ~u m mv t -v,„) 

i,3=i n>k>m>l 

x(he(h)r(p^_ in _ m (f(u n ) ■ ..f( Ul )),e(v n ) ■ ..e( Vl )). (46) 

Using the expression for the projection of the product of the total currents (|45p and 
Formula (1281) we obtain 



Z in) (u n , . . . ,ui,v n , . . .,vx] A) 



<t(!)>ct(!') 

fc>m ~ *V( m ))0(-A - (m - l)ft) 

TT 6>(^fc - V m - ft) ST TT °(Vg(l) ~ V<r{V) + ^) 

<t(0><t(!') 



X 



n a( L k\Tf at \ TT ^ m ~ V«0 ~ * - (m - l)ft)O(ft) 
0{u k - v a{m) + ft)[[ 6(u k - v a(m) ) H at_ x _t m _ 1)h) 

k>m k<m m=l 



From this formula we see that the expression (j4"7j) defines a holomorphic function of the 
variables u%. 

Theorem 1 The set of functions {Z^ n \u n , . . . , Ui, v n , . . . , Vx\ A)} n>1 defined by For- 
mula Ub]) satisfies the conditions of Propositions^ (H and the initial condition 
They coincide with the partition functions of the SOS model with DWBC: 

Z+~(u n , . . . ,ui, v n , . . .,vi, A) = 

nO(v k -Vm-ft)^ yj 9(v a{l) - V a{ i>) + ft) TT Q( v 
Tt r~ 1^ II Tt fcT 11 ¥*-"»W 

n>k>m>l ff£6„ i<i' v ' v ' l<fe<m<n 

<r(i)>CT(!') 
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TT n( , *-\ TT A(Mm ~ ^ ~ A - (m - l)fi)g(g) 

n>fc>m>l m=l v v 11 

The initial condition (I22p is satisfied because Formula (I34p is satisfied. The first factor 
in the right-hand side of (j46p is symmetric with respect to both sets of variables. The 
symmetry with respect to the variables {u} and the variables {v} follows from the com- 
mutation relations (|26|) and (1251) respectively. Formula (147]) implies that (146|) are elliptic 
polynomials of degree n with character ( fl8l) in the variables Mj, in particular in u n . We 
now substitute u n = v n — h to (j47j). The non- vanishing terms in the right-hand side 
correspond to the permutations a G S n satisfying a(n) = n. Substituting u n = v n — h~ 
into these terms, one obtains the recursion relation (TTTJT) . □ 



5. Degeneration of the partition function 



In this section, we investigate the trigonometric degenerations of the formulae obtained 
in the elliptic case. In particular, taking the corresponding trigonometric limit in the 
expression for the SOS model partition function (|48j) reproduces the expression for the 
6- vertex partition function ([5]). 

First we consider the degeneration of the i?-matrix, the matrix of Boltzmann 
weights, which defines the model. To do so we need the formula for the trigonometric 
degeneration (r — > ioo) of the odd theta function, 

sin 7tu 



lim 9{u) 



Ti 



In terms of the multiplicative variables z 
as follows: 



27rie ni{u+v) lim 0( 



u 



e 27rm w 



W . 



this formula can be rewritten 



Multiplying the i?-matrix ffTU]) by 2nie m ( u+vS) and taking the limit we obtain the following 
matrix which depends rationally on the multiplicative variables z, w and on the 
multiplicative parameters q = e mh , /i = e 2nlX : 

R(z, w; n) = 2irie ni{u+v) lim R(u 



( zq — wq~ 



\ 











(z-w)(fiq-q- 1 ) 





v;X) 




(z-w^)(q-q 1 ) 
(z-w)(nq 1 -q) 











(49) 



(zq — wq ) / 

The matrix (14^1) inherits the property of satisfying the dynamical Yang-Baxter equation 
and it defines a statistical model called the trigonometric SOS model. 



R(z, w) = lim R(z,w;fi) 



.(50) 
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To obtain the non-dynamical trigonometric case we need to implement the 
additional limit A — > — ioo implying fi — > oo (or A — > ioo implying fi — > 0): 

/ zg-wg- 1 \ 

g(z — tu) (g — g _1 )u; 

(q — q~ x )z q^ 1 (z — w) 

V zq - wq^ 1 J 

The matrix (!50j) differs from the matrix of the Boltzmann weights of the 6-vertex 
model (j2J) by the transformation (12 Taking into account Remark H] (from Subsection 
3.2), we conclude that both matrices ([2]) and (1501) define the same partition function 
Z({z},{w}) with DWBC0. 

To obtain the partition function with DWBC for the trigonometric SOS model, 
one should multiply the partition function with DWBC for the elliptic SOS model by a 
factor and take the trigonometric limit: 

■n 

Z+-({z},{w} ]f i)= TT (2rcie^ u ^) lim Z+-({u},{v};\) = 

k,j=l 

n>k>m>X a&on Ki' v ' v ' 

<r(i)>o-(!') 

(z m - «; < T( m )Aig 2(ro " 1) )(g - g" 1 ) 



X 



I [ {z k q-w a{m) q l ) j [ (z fc - w CT ( m) ) 



(1 - aq 2 ^- 1 )) 

n>k>m>l Kk<m<n m=l v y 



It is easy to prove that Formula (jSJ) is obtained from Formula (15 ip by taking the limit 
Z({2},{w})= lim Z^({4,{ W };//). 
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Appendix A. Interpolation formula for elliptic polynomials 

A group homomorphism T — > C x , where T = Z + rZ and C x is the multiplicative 
group of nonzero complex numbers is called a character. Each character x an d integer 
number n define a space O n (x) consisting of the holomorphic functions on C with the 
translation properties 

0( M + 1) = <t>{u + r)= X (r)e- 2 ™™>H. 

If n > then dim0 n (x) = n (and dimO n (x) = if n < 0). The elements of the space 
n (x) are called elliptic polynomials (or theta-functions) of degree n with character \- 

Proposition 4 Let {0,}™ =1 be a basis ofQ n (x), with character x(l) = (— l) n , x( T ) = 
= (— l) n e 2ma , then the determinant of the matrix \\<f>j(ui)\\<i,j<n is equal to 

n 

det H^tOH = C ■ 0(J> fe - a)l[6(u t - Uj ), (A.l) 

k=l i<j 

where C is a nonzero constant. 
Consider the ratio 

det 110,(^)11 



0(ELi^-a)n^K-%)' (A,2) 

This is an elliptic function of each Ui with only simple poles in any fundamental domain 
(the points Ui satisfying Y^l=\ u k — & £ T). Therefore, it is a constant function of each -Uj. 
Thus this ratio does not depend on Ui and we have to prove that it does not vanish, that 
is that the determinant det ||0j(-Uj)|| is not identically zero. Let us denote by A**''"'**: 
the minor of this determinant corresponding to the ii-th, . . ., ifc-th rows and the Ji-th, 
. . ., jfc-th columns. Suppose that this determinant is identically zero and consider the 
following decomposition 

n 

det ||^)|| = ^(-i) fe+ V,(yi)At::;l l _ lifc+1 ,..., n . (a.3) 



k=l 



Since the functions <pk{yi) are linearly independent, the minors A 2 ''"'^ k+1 n are 



identically zero. Decomposing the minor A^'"'™ we conclude that the minors 



^2, 'S-i,k+i,...,n are identically zero, and so on. Finally, we obtain that A™ = 4> n {yn) 
is identically zero which cannot be true. □ 

Lemma 2 Let us consider two elliptic polynomials Pi, P 2 G Q n {x)> wherexi^) = (— l) n ; 
x( r ) — (— l) n e a , and n points Ui, i = l,...,n, such that Ui — Uj T , i ^ j, 
and Y^k=i u k ~ a & r. If the values of these polynomials coincide at these points, 
Pi{ u i) — ^2(^1), then these polynomials coincide: P\{u) = Pziu). 

Decomposing the polynomials under consideration as P a {u) = Yl^iPa^ii 11 ), a = 1>2, 
we obtain the system of equations 

n 

1=1 
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with respect to the variables p\ 2 = p\ — p\. We have proved that the determinant of 
this system is equal to flA.lj) and therefore is not zero. Hence, this system has only the 
trivial solution p\ 2 = 0, but this implies P\{u) = P2{u). □ 
Let P e On(x) be an elliptic polynomial, where = ( — 1)™; x( T ) — (—l) n e 2ma , 
and Ui, i — 1, . . . , n, be n points such that U{ — Uj T , i ^ j, and Y^k=i u k — a T. 
This polynomial can be recovered from the values at these points: 

P(u) = ± PM ?!* - u + a ~ S=| <1 JJ pZ*L (A.4) 

Indeed, the right hand side belongs to Q n (x), this equality holds at the points u = U{. 
Using Lemma [21 we conclude that (IA.4I) holds at all u £ C. 
Consider the meromorphic functions 

M _ 9{u,j -u + X- (n-2j + 2)h) 6{u k - u - h) 
Qj {U} ~ 9( Uj -u + h) \} 2 6(u k -u + h)' 

It is easy to check that the functions 

n 

p A u ) = n °( uk ~ u+ n )QA u ) = 

k=2 

j-1 n 

= 6(uj - u + X-(n — 2j + 2) ft) JJ 9{u k -u-h) JJ 6(u k -u + H) 

k=2 k=j+l 

belong to e n _!(x), where = (-l) n_1 , x(^) = (-l)™"^ 2 ™, a = A + EJU«*- 
Since A T, the polynomials -Pj(m) can be recovered from by their values Pj(ui) via the 
interpolation formula (1A.4I) . Taking into account the relation between Qj{u) and Pj(u) 
we obtain Formula fl38l) . □ 



i?e/erences 

[1] Khoroshkin S and Pakuliak S 2005 Theor. and Math. Phys. 145 1373-99 

[2] Enriquez B and Felder G 1998 Commun. Math. Phys. 195 651-89 

[3] Enriquez B and Rubtsov R 1997 Ann. Sci. Ecole Norm. Sup. 30 821-46 

[4] Enriquez B and Rubtsov V 1999 Israel J. of Math. 112 61-108 

[5] Enriquez B and Rubtsov V 1999 Amer. Math. Soc. Transl. Ser. 191 33-65 

[6] Drinfeld V 1988 Soviet Math. Doc. 36 212-16 

[7] Ding J, Khroshkin S and Pakuliak S 2000 Theor. and Math. Phys. 124 1007-37 

[8] Enriquez B, Khoroshkin S and Pakuliak S 2007 Commun. Math. Phys. 276 691-725 

[9] Izergin A 1987 Soviet Phys. Dokl. 32 878-79 

[10] Rosengren H 2008 An Izergin-Korepin-type identity for the 8VSOS model, with applications to 
alternating sign matrices Preprint arXiv:0801.1229 [math. CO] 

[11] Tarasov V and Varchenko A 1997 Asterisque 246 

[12] Khoroshkin S, Pakuliak S and Tarasov V 2007 J of Geom. and Phys 57 1713-32 
[13] Felder G and Schorr A 1999 J. Phys. A: Math. Gen. 32 8001-22 

[14] Date E, Jimbo M, Kuniba A, Miwa T and Okado M (1987) Nuclear Phys. B 290 231-73 

* We can require the condition m — Uj $ T, because the Mi's in Formula (|38p are formal variables. 



SOS model partition function and the elliptic weight functions 



21 



[15] Pakuliak S, Rubtsov V and Silantyev A 2007 Classical elliptic current algebras Preprint 
arXiv:0709.3592 [math.QA] Accepted for publication in J. of Generalized Lie Theor. and App. 
[16] Kuperberg G 1996 Internat. Math. Res. Notices 10 139-50 
[17] Odesskii A and Feigin B 1997 Internat. Math. Res. Notices 11 531-39 

[18] Felder G, Varchenko A and Tarasov V 1997 Topics in singularity theory Amer. Math. Soc. Transl. 

Ser. 180 45-75 (Providence: AMS) 
[19] Enriquez B, Pakuliak S and Rubtsov V Contemporary Math. 433 177-90 



